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ABSTRACT 

We present a new determination of the cluster mass function and velocity dispersion function in a 
volume ^\Q^h^y[^c~^ using data from the Fourth Data Release of the Sloan Digital Sky Survey (SDSS) 
to determine virial masses. We use the caustic technique to remove foreground and background galaxies. 
The cluster virial mass function agrees well with recent estimates from both X-ray observations and 
cluster richnesses. Our determination of the mass function lies between those predicted by the First- 
Year and Three- Year WMAP data. We constrain the cosmological parameters and erg and find good 
agreement with WMAP and constraints from other techniques. With the CIRS mass function alone, we 
estimate = 0.24t°;J^ and erg = 0.92j:|5;^^, or cts = 0.84±0.03 when holding = 0.3 fixed. We also 
use the WMAP parameters as priors and constrain velocity segregation in clusters. Using the First and 
Third- Year results, we infer velocity segregation oiugxy/crDM ~0.94±0.05 or 1.28±0.06 respectively. The 
good agreement of various estimates of the cluster mass function shows that it is a useful independent 
constraint on estimates of cosmological parameters. We compare the velocity dispersion function of 
clusters to that of early-type galaxies and conclude that clusters comprise the high-velocity end of the 
velocity dispersion function of dark matter haloes. Future studies of galaxy groups are needed to study 
the transition between dark matter haloes containing individiual galaxies and those containing systems 
of galaxies. The evolution of cluster abundances provides constraints on dark energy models; the mass 
function presented here offers an important low redshift calibration benchmark. 

Subject headings: galaxies: clusters — galaxies: kinematics and dynamics — cosmology; observations 



1. INTRODUCTION 

Clusters of galaxies are the most massive gravitationally 
relaxed systems in the universe, so the observed cluster 
mass function is a sensitive probe of cosmological param- 
eters. Galaxy cluster abundances are most sensitive to 
the matter density of the universe and cg, the rms 
fluctuations in spheres of radius 8/i~^Mpc and the usual 
normalization of the linear power spectrum (Henry & Ar- 
naud 1991; BahcaU & Gen 1993). The evolution of the 
mass function is an excellent probe of the growth of struc- 
ture, and thus can constrain the properties of dark en- 
ergy (e.g., Haiman et al. 2001; Hu 2003; Vikhlinin et al. 
2003; Henry 2004; Majumdar & Mohr 2004; Voevodkin & 
Vikhlinin 2004). Robust constraints from cluster evolu- 
tion are only possible with an accurate determination of 
the cluster mass function in the nearby universe. 

Most recent estimates of erg from the local cluster mass 
function find values of crg«0.6-0.8 for 57^=0.3 (Pierpaoli 
et al. 2001; Reiprich & Bohringer 2002; Seljak 2002; 
Viana et al. 2002; Pierpaoli et al. 2003; Schuecker et al. 
2003; Viana et al. 2003; BahcaU et al. 2003; Voevodkin 
& Vikhhnin 2004; Eke et al. 2005; Dahle 2006) (there 
is a strong degeneracy between erg and 17™). Estimates 
based on cosmic shear usually find larger values of (TgR:i0.8- 
1.0 (Massey et al. 2005; Hoekstra et al. 2005; Semboloni 
et al. 2006), although some investigators find smaller val- 
ues (Hamana et al. 2003; Heymans et al. 2005). Re- 
cent attempts to model the relation between dark mat- 



ter haloes and galaxies (van den Bosch et al. 2003; Tinker 
et al. 2005) suggest that smaller values of CTg«0.6-0.8 are 
required to match cluster mass-to-light ratios (Carlberg 
et al. 1996; Bahcall et al. 2000; Rines et al. 2004). Various 
estimates based on bulk flow models seem to be converg- 
ing on intermediate values (see Table 3 of Pike & Hud- 
son 2005). The WMAP satellite has produced tight con- 
straints on many cosmological parameters (Spergel et al. 
2003, 2006). Among the largest changes between the First- 
Year and Three- Year WMAP results (hereafter WMAPl 
and WMAP3) was a decrease in both erg and Q,m (Spergel 
et al. 2006). The WMAP3 results provide a much bet- 
ter fit to the mass function of X-ray clusters from the 
ROSAT Highest X-ray FLUx Galaxy Cluster Survey (HI- 
FLUGCS, Reiprich & Bohringer 2002) than the WMAPl 
results (Reiprich 2006). 

The greatest difficulty in measuring the cluster mass 
function is obtaining sufficiently accurate mass estimates. 
Three well-established methods exist for measuring clus- 
ter masses. The dynamics of cluster galaxies provided the 
first evidence for dark matter when Zwicky (1933, 1937) 
applied the virial theorem to the Coma cluster. The sec- 
ond method is to use the properties of the hot intracluster 
medium (ICM), whose distribution and temperature probe 
the gravitational potential of the cluster. Finally, gravi- 
tational lensing can provide very accurate mass estimates 
in the centers of clusters and somewhat less accurate esti- 
mates at large radii. 

All three techniques have been used extensively, but 
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they are each subject to potentially large systematic un- 
certainties. There may be velocity bias between cluster 
galaxies and dark matter particles, although the magni- 
tude and direction of this bias is under debate (e.g., KaufF- 
mann et al. 1999b; Dicmand ct al. 2004; Faltcnbachcr 
et al. 2005; Benatov et al. 2006). Chandra and XMM- 
Newton observations show that the physical properties of 
the ICM are complex and involve significant interaction 
with AGN (e.g., McNamara et al. 2000; Markevitch et al. 
2000; Fabian ct al. 2000; VikhUnin ct al. 2001; Fabian ct al. 
2006). Finally, mass estimates from gravitational lensing 
(especially at large radii) are subject to significant uncer- 
tainties due to lensing by other structures along the line- 
of-sight (e.g., Metzler et al. 2001), although these might be 
overcome by combining information from strong and weak 
lensing (Natarajan et al. 2002; Bradac et al. 2004). These 
potential systematics have led some investigators to invert 
the problem and determine cluster properties by match- 
ing the mass function predicted by a cosmological model 
to the observed richness function or luminosity function 
(Berlind et al. 2006; Stanek et al. 2006). 

Despite these challenges, the cluster mass function is the 
subject of much study. Most recent measurements of the 
cluster mass function have been made with either X-ray 
data (Reiprich & Bohringer 2002; Allen et al. 2003; Vo- 
evodkin & Vikhlinin 2004) or optical richness data (Bah- 
call et al. 2003). The most significant uncertainty in deter- 
mining (Tg from X-ray observations is the normalization of 
scaling relations, either the M—Tx or the M—Lx relation. 
The normalizations are often determined from hydrody- 
namical simulations (e.g., Borgani et al. 2001), but differ- 
ent simulations produce a range of normalizations (Pier- 
paoli et al. 2003; Viana et al. 2003). Another potential 
concern in determining scaling relations is Malmquist bias 
(Stanek et al. 2006). Voevodkin fc Vikhlinin (2004) avoid 
these scaling relations by estimating the mass of the intr- 
acluster medium (ICM) and measuring the baryonic mass 
function (a similar analysis was presented by Shimasaku 
1997). This method requires assumptions about: the rela- 
tive contribution of stars and gas to the total baryon mass, 
the ratio of the baryon fraction in clusters to the global 
value, and the mass dependence of this ratio. A potentially 
large systematic uncertainty in all of these X-ray studies is 
a possible offset between the temperature Tspec measured 
by X-ray satellites and the emission- weighted temperature 
Teu, calculated in hydrodynamical simulations (Rasia ct al. 
2005). This effect can be understood as the excess contri- 
bution of line emission from cooler gas in an ICM with a 
variety of temperatures, thus causing Tspec to be an under- 
estimate of Tew Correcting for this effect would increase 
the estimated cluster masses and thus increase the best-fit 
values of erg (Rasia et al. 2005; Shimizu et al. 2006). A sim- 
ilar systematic effect is suggested by the higher normaliza- 
tion of the mass-tcmpcrature relation found by Vikhlinin 
et al. (2006) from a study of relaxed clusters with temper- 
ature profiles from Chandra data. An increase in this nor- 
malization would increase the inferred cluster masses and 
hence the overall mass function, shifting the constraints 
towards higher values of and (Ts. 

Viana et al. (2002) use a hybrid approach of combin- 
ing the X-ray luminosity function of clusters with a mass- 
luminosity relation determined from weak lensing (Sheldon 



et al. 2001). They find results consistent with other X-ray 
studies, but note that large systematic uncertainties are 
possible. Recently, Dahle (2006) used weak lensing mass 
estimates of an X-ray selected cluster sample to measure 
the mass function. He found good agreement with many 
of the constraints from X-ray estimates but lower values 
of Om and erg than recent cosmic shear studies. 

To our knowledge, the most recent cluster virial mass 
function is that of Girardi et al. (1998a). The advent of 
large-scale; rcxlsliift surveys, especially the Sloan Digital 
Sky Survey (SDSS, Stoughton et al. 2002), enables more 
accurate virial mass estimates due to their size and unifor- 
mity. Compared to X-ray studies, virial mass estimates are 
sensitive to larger scales (r2oo rather than rsoo, where ta 
is the radius within which the enclosed density is A times 
the critical density), which allow for comparison with the- 
oretical mass functions with significantly less extrapola- 
tion (White 2002). Also, virial masses can be estimated 
for poor clusters and rich groups, whereas X-ray mass es- 
timates of these systems are complicated by possible en- 
ergy input from supernovae and AGN (e.g., Loewenstein 

2000) . Probing these smaller masses enables a direct con- 
straint on fluctuations on the scale Sft.^'^Mpc, rather than 
the ~14/i~^Mpc scales probed by ~lO^^/i~^M0 clusters 
(Pierpaoli et al. 2001). Eke et al. (2005) used a simphfied 
version of the virial theorem to estimate the group mass 
function from an optically-selected group catalog in the 
Two Degree Field Galaxy Reclshift Survey (CoUess et al. 

2001) . Differences in our analysis include: X-ray selec- 
tion rather than optical selection of clusters, use of the 
full virial theorem with corrections for the surface pres- 
sure term, much better sampling of individual systems 
(the average number of cluster members is 50), and more 
conservative rejection of interlopers (thus significantly re- 
ducing scatter in the mass estimates). Estimates of the 
group mass function using virial masses are usually lim- 
ited by systematics in groiip selection, mass estimation, 
and cosmic variance (Girardi & Giuricin 2000; Martmez 
et al. 2002; Heinamaki et al. 2003). 

Here we use the virial theorem to estimate cluster masses 
of an X-ray selected sample of clusters with data from the 
Sloan Digital Sky Survey. The use of X-ray selection re- 
duces the impact of projection effects on cluster selection. 
Perhaps the greatest advantage of X-ray selection of our 
sample is that the selection function is well understood 
and can be computed directly. The cluster sample we use 
is the Cluster Infall Regions in SDSS (CIRS) sample of 
Rines & Diaferio (2006) with some minor modifications 
described below. 

We describe the data and the cluster sample in § 2. 
In § 3, we estimate the mass function using both virial 
masses and caustic masses and discuss potential system- 
atic effects. We compute the velocity dispersion function 
of clusters in §4 and compare it to the velocity dispersion 
function of early-type galaxies in SDSS. We discuss our re- 
sults and conclude in §5. We assume Ht^ = lOO/i km s~"^, 
and a fiat ACDM cosmology (Oa = 1 — ^m) throughout. 
Where not stated explicitly, we assume VLm — 0.3 for initial 
calculations. 

2. THE CIRS CLUSTER SAMPLE 

2.1. Sloan Digital Sky Survey 
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The Sloan Digital Sky Survey (SDSS, Stoughton et al. 
2002) is a wide-area photometric and spectroscopic sur- 
vey at high Galactic latitudes. The Fourth Data Release 
(DR4) of SDSS includes 6670 square degrees of imag- 
ing data and 4783 square degrees of spectroscopic data 
(Adelman-McCarthy et al. 2006). 

The spectroscopic limit of the main galaxy sample of 
SDSS is r=17.77 after correcting for Galactic extinction 
(Strauss et al. 2002). Assuming the luminosity function of 
Blanton ct al. (2003), the spectroscopic limit of SDSS cor- 
responds to Mo*.i^-hl at z = 0.092. The Cluster And Infall 
Region Nearby Survey (CAIRNS Rines et al. 2003) found 
that well-sampled clusters contain dense envelopes known 
as caustics (Diafcrio 1999) which provide a direct means 
of determining cluster membership. CAIRNS sampled all 
galaxies brighter than about M* + 1 and often fainter. Fig- 
ure 7 of Rines et al. (2005) shows the infall patterns for 
galaxies brighter than M'^ -f 1 in the CAIRNS clusters. 
These patterns are readily apparent, but less well-defined 
than the full CAIRNS samples (Rines et al. 2003). Thus, 
we expect that infall patterns of clusters with masses sim- 
ilar to CAIRNS clusters should be apparent to z <0.1, 
though not much further. Rines & Diaferio (2006) show 
the infall patterns of the clusters at z <0.1 in SDSS and 
compute the mass profiles from the caustic technique and 
the virial theorem. 

Note that the SDSS Main Galaxy Survey is ~85-90% 
complete to the spectroscopic limit. The sinvcy has w7% 
incompleteness due to fiber collisions (Strauss et al. 2002), 
which are likely more common in dense cluster fields. Be- 
cause the target selected in a fiber collision is determined 
randomly, this incompleteness can theoretically be cor- 
rected for in later analysis. From a comparison of SDSS 
with the Millennium Galaxy Catalogue, Cross et al. (2004) 
conclude that there is an additional incompleteness of ~7% 
due to galaxies misclassificd as stars or otherwise missed 
by the SDSS photometric pipeline. For our purposes, the 
incompleteness is not important provided sufficient num- 
bers of cluster galaxies do have spectra. 

2.2. X-ray Cluster Surveys 

Because SDSS surveys primarily low-redshift galaxies, 
the best sampled clusters are both nearby and massive. 
We therefore search X-ray cluster catalogs derived from 
the ROSAT All-Sky Survey for clusters in DR4. RASS 
(Voges ct al. 1999) is a shallow survey but it is sufhciently 
deep to include nearby, massive clusters. RASS covers vir- 
tually the entire sky and is thus the most complete X-ray 
cluster survey for nearby clusters. 

Published cluster catalogs derived from the RASS in- 
clude the X-ray Brightest Abell Cluster Survey (XBACS 
Ebeling et al. 1996), the Bright Cluster Survey and its ex- 
tension (BCS and eBCS Ebehng et al. 1998, 2000), the 
Northern ROSAT All-Sky galaxy cluster survey (NORAS 
Bohringer ct al. 2000), and the ROSAT-ESO flux limited 
X-ray galaxy cluster survey (REFLEX Bohringer et al. 
2001). We refer the reader to the individual catalog pa- 
pers for detailed descriptions of the construction of the 
catalogs (Rines & Diaferio 2006, contains a summary). 
When combined, these catalogs cover virtually the entire 
sky at high Galactic latitudes {\b\ >20°) to a flux limit of 
Ri3xl0~^^erg cm~^ s~^. None of these catalogs is com- 



plete, so by combining them, we create a substantially 
more complete composite catalog. 

2.3. Definition of the Mass Function Sample 

When multiple X-ray fluxes are available for a cluster, 
we use the most recently published value. The order of 
preferences is therefore: REFLEX, NORAS, BCS/eBCS, 
XBACs. The details of the flux determinations in the dif- 
ferent catalogs vary somewhat: NORAS and REFLEX 
both measure fluxes with Growth Curve Analysis, while 
BCS/eBCS uses Voronoi Tessellation and Percolation. De- 
spite these differences, the derived fluxes agree fairly well 
(see Figure 21 of Bohringer et al. 2000, for a direct compar- 
ison of NORAS and BCS/eBCS), so combining the cata- 
logs is a reasonable procedure. By using the most recently 
published flux value, we use the GCA flux if available and 
the VTP flux only when no GCA flux is available. 

We define a flux-limited and redshift-limited sample of 
clusters with the criteria /x > 3 x 10~^^erg s^^cm"^ (0.1- 
2.4 keV) and z <0.10. Note that variations in the method 
of determining flux in different catalogs may affect the 
precise flux limit. The only cluster in XBACS but not in 
any other catalog is Al750b: XBACS identifies two X-ray 
sources (A1750a/b), whereas other catalogs treat them as 
one source. We follow Bohringer et al. (2004) and treat 
A1750 as a single source rather than two separate sources 
as in XBACs. Similarly, the galaxy NGC5813 is boimd to 
the NGC5846 group (Mahdavi et al. 2005). Because the 
dynamics of NGC5813 are dominated by the NGC5846 
system, we eliminate NGC5813 from the sample (both are 
eliminated by our minimum redshift cutoff below). We 
inspect the redshift data aroimd each cluster to confirm 
the cluster redshift and find that A2064 has an incorrect 
redshift (and X-ray luminosity) listed in NORAS. The cor- 
rect redshift is 0.0738 instead of 0.1076. We correct the 
X-ray luminosity accordingly. Our final flux and volume 
limited sample contains 72 clusters within the SDSS DR4 
spectroscopic footprint. We will refer to this sample as the 
CIRS (Cluster Infall Regions in SDSS) clusters in the rest 
of this paper. The completeness of the CIRS sample is 
limited by the completeness of the underlying cluster cat- 
alogs. However, by combining clusters from the various 
catalogs we should be more complete than any one of the 
catalogs. The local nature of the sample {z <0.1) likely 
improves the completeness of the CIRS sample relative 
to a flux-limited sample. We discuss the modest poten- 
tial incompleteness in more detail in §3.3. The clusters 
are an unbiased sample: the selection of the CIRS sam- 
ple is based purely on X-ray flux and the footprint of the 
SDSS DR4 spectroscopic survey. We confirm this claim 
with a V/Vmax test: we find < V/Vmax >=0.518±0.035 
compared to an expected value of 0.5 for a complete, uni- 
form sample (see §3.3). A homogeneous, complete, all-sky 
catalog of extended X-ray sources in RASS would allow a 
significant improvement in the determination of the mass 
function. Unfortunately, no such catalog currently exists. 

For computing the mass function, we omit three clusters 
(Virgo, NGC4636, and NGC5846) with z < 0.02 because 
their peculiar velocities are likely large compared to their 
redshifts. From our analysis of individual systems in CIRS 
(see §6 of Rines & Diaferio 2006), we eliminate A1035A/B 
and A1291A/B as superpositions which would lie below 
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the flux limit if they were resolved in the catalogs. Fi- 
nally, we eliminate A2249 because its center lies ~2' off 
the edge of the SDSS DR4 spectroscopic footprint. After 
eliminating these six clusters, our mass function sample 
includes 66 clusters. We test the effects of returning these 
clusters to the mass function sample in §3.3. 

3. THE CLUSTER MASS FUNCTION 

3.1. Estimating the Mass Function 

We compute both the standard cluster mass function 
dn{M) / dlogM and the cumulative mass function n(>M), 
the number density of clusters more massive than M. To 
compute the mass functions, we use the 1/Vmax estimator 
(Schmidt 1968), where Vmax{Lx) is the maximum (comov- 
ing) volume a cluster with X-ray luminosity Lx would lie 
within our flux and redshift limited sample. In each loga- 
rithmic mass bin, we sum the clusters 

dn{M) _ 1 ^ 1 

dlogM dlogM ^ Vma4Lx,i) ^ ' 

where the sum is over clusters within the mass bin. By 
using Vmax{Lx) rather than Vmax{M), we avoid the need 
to know the slope, normalization, and scatter of the scal- 
ing relation Lx — M to properly calculate Vmax (Reiprich 
& Bohringer 2002). Figure 1 shows the maximum volume 
probed by CIRS as a function of X-ray luminosity. Al- 
though CIRS covers a much smaller area than HIFLUGCS, 
Vmax{Lx) is larger for CIRS for \ogLx < 43.8 because 
of the smaller flux limit. We calculate Vmax assuming a 
flat n„i=0.3 cosmology. We repeated our analysis using 
f^m=0.1 and Q,fn=^ to calculate Vmax and find that this 
has a negligible effect due to the local nature of the sample. 
We correct our mass functions for a 5% incompleteness (see 
§3.3.1) which we assume applies uniformly to the sample. 

We estimate the uncertainty in the mass function as 
[A((in(M)/dlogM)]2 = J:^[^/Vmax{Lx,^n Figurc 2 
shows the CIRS mass function and Figure 5 shows the 
cumulative mass function. We compute the masses in 
two ways, one using the virial theorem (red lines in Fig- 
ures 2 and 5) and the other (blue lines) using the caustic 
technique (Diaferio 1999; Rines & Diaferio 2006). We use 
the caustic technique to identify cluster members in both 
techniques. As a membership classification algorithm, the 
caustic tcchniqiic is similar to the "shifting gappcr" tech- 
nique (e.g., Fadda et al. 1996), but is less sensitive to in- 
terlopers. The mass functions computed with these two 
mass estimators agree quite well, a confirmation that the 
caustic technique provides reliable mass estimates (Diafe- 
rio 1999: Gcllcr ct al. 1999; Rines ct al. 2000; Rcisenegger 
et al. 2000; Rines et al. 2001, 2002; Biviano & Girardi 
2003; Rines et al. 2003; Rines & Diaferio 2006). We esti- 
mate the masses M200 within r200) where ta is the radius 
within which the average density is Apc, where pc is the 
critical density. We determine r2oo directly from the non- 
parametric mass profiles. 

For the virial masses, we obtain an initial estimate of 
?'200,t)ir from the virial mass profiles (Rines & Diaferio 
2006). We then apply a correction factor of 8% to ac- 
count for the surface pressure term of the virial theorem. 
This factor is calculated from Equation 8 of Girardi et al. 
(1998b) assuming galaxies are on isotropic orbits and an 



NFW mass profile (Navarro et al. 1997) with a concentra- 
tion parameter of c = r2oo /^s = 5 (r^ is a scale radius in 
the NFW profile). The assumption of isotropic orbits is 
supported by numerous observations (van der Marel et al. 
2000; Rines ct al. 2003; Lokas & Mamon 2003; Biviano & 
Katgert 2004; Rines & Diaferio 2006; Benatov et al. 2006; 
Lokas et al. 2006), as is the value of C200 = 5 (Markevitch 
et al. 1999; Lin et al. 2004; Rines & Diaferio 2006). After 
this correction, the enclosed density within r2on,vir is now 
178pc. We therefore reduce the mass by a further 3.3% to 
obtain our final estimate of M2oo,vir- The concentration 
c is expected to be a weak function of mass, c oc M~°-^^ 
(Bullock et al. 2001). This corresponds to a 25-30% change 
over the mass range we probe. This change has little effect 
on the mass estimates; assuming c=5 leads to a mass un- 
derestimate (overestimate) of 3% (5%) for a cluster with 
c =10 (2.5). The CIRS clusters are quite well sampled, 
with an average of ss50 members per cluster. The statisti- 
cal uncertainties in M^r are on average 13%. The caustic 
and virial mass functions are in good agreement, further 
confirming the good agreement between caustic and virial 
mass estimates (Rines et al. 2003; Diaferio et al. 2005; 
Rines k Diaferio 2006). 

The original prescription for calculating the mass func- 
tion was proposed via the Press-Schechter formalism 
(Press & Schechter 1974). Numerical simulations pre- 
dict relatively more massive systems and fewer less mas- 
sive systems than Prcss-Schcchtcr theory (Slicth & Tor- 
men 1999). Jenkins et al. (2001) provided fitting formu- 
lae for a universal mass function that can be evaluated 
for many cosmological models. In particular, their mass 
function accurately reproduces the mass function of dark 
matter haloes in the Hubble Volume simulation. Sev- 
eral investigators conclude that the mass function pro- 
vided by Equation (B3) of Jenkins et al. (2001) is close 
to a universal mass function and can thus be used to con- 
strain cosmological parameters (White 2002; Evrard et al. 
2002; Hu & Kravtsov 2003). For this formula, haloes are 
placed at the most bound particles, then the radius of the 
halo is increased until the enclosed spherical overdensity 
is 180 times the background (not critical) density. This 
mass M18O6 (the "b" in the subscript indicates background 
rather than critical density) must then be converted to 
M200 assuming an NFW profile with c=5 (White 2002; 
Evrard et al. 2002; Hu & Kravtsov 2003). Note that X- 
ray mass estimates often require conversion to M500; the 
conversion to M200 is less of an extrapolation and less sen- 
sitive to the assumed concentration parameter. 

Figure 2 shows the mass functions for the best-fit cos- 
mological parameters from WMAPl (Spergel et al. 2003) 
and WMAP3 (Spergel et al. 2006) at z=0.062 (the mean 
redshift of our sample). The WMAP mass functions in 
Figure 2 are convolved with our mean mass uncertainty 
'^logM ~ 0.056 for appropriate comparison with the ob- 
served mass functions. Some of the largest differences be- 
tween the WMAPl and WMAP3 were in ag and Ctm, the 
two parameters which have the largest impacts on the mass 
function. Figure 2 shows that the CIRS mass function 
lies between the WMAPl and WMAP 3 results. Reiprich 
(2006) showed that the HIFLUGCS mass function esti- 
mated using the X-ray properties of clusters lies closer 
to the WMAP3 results (Reiprich & Bohringer 2002), al- 
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though other mass functions derived with X-ray mass trac- 
ers find values closer to the WMAPl results (Pierpaoli 
et al. 2003; Viana et al. 2003). 

The dotted lines in Figure 2 show the mass func- 
tion computed without imposing a minimum rcdshift 
on the sample. Clearly, the systematic uncertainty 
in the mass function due to cosmic variance becomes 
large at M<10^'^h~^MQ (see also Figure 1). We mini- 
mize for the CIRS mass function in the mass range 
logA'/2oo=[13.9,15.1] by calculating the Jenkins ct al. 
(2001) mass function for a given pair of {flmjCrs) and shift- 
ing the mass scale from Misob to M2oo- We set F according 
to 

F(0„,/i) = fi„/if^) e^p{-nt,-^/2h^) (2) 

(Sugiyama 1995) with Tq = 2.726ii:, h = 0.7 and Qf, = 
0.0223/1^ (Spergel et al. 2006). Holding r=0.21 fixed does 

not strongly affect our results. We convolve the mass func- 
tion with our mean mass uncertainty trj^g^^^ = 0.056 ac- 
cording to 



dn{M)/dM 



x(27raf,,M)"'/' 



exp 



1 r+°° 



dM 

(logM' - logM)2 

■^"logM 



d\ogM{3) 



where Vmax (M) is computed using the scaling relation 

log(M2oo/i/M0) = 0.7631og(Lx,44) + 14.62 (4) 

from Rincs & Diaferio (2006) where Lxa4 is the X-ray 
luminosity in units of 10^^/i^^crgs^^. This relation has a 
different normalization than the relation found by Popesso 
et al. (2005), but is consistent with the relation of Reiprich 
& Bohringer (2002). The weighting by Vmax is required 
because the CIRS mass function is flux- limited for z < 0.1. 
We repeat the fits without weighting by Vmax and find neg- 
ligible differences. At smaller masses, the volume probed 
is fairly small (Figure 1), so the mass function estimate is 
subject to significant uncertainty due to cosmic variance 
(Hu & Kravtsov 2003). In addition, the number of clus- 
ters at these low masses is probably not large enough to 
sample the scatter in Lx- Figure 3 shows the 68% and 
99.7% confidence levels for 0„i and inferred from the 
virial mass function (solid contours) and the caustic mass 
function (dashed contours). With the CIRS virial mass 
function, we estimate flm = 0.24lg;j!,g and erg = 0.92+o;^g, 
or fT8 = 0.84±0.03 when holding iim = 0.3 fixed (all as- 
suming h=0.7). The open squares in Figure 3 show the 
WMAPl and WMAP3 estimates with their 68% uncer- 
tainties as crrorbars (Spergel et al. 2006, Table 2). The 
CIRS mass function is consistent with the WMAP3 results 
within their 95% confidence levels, while the WMAPl re- 
sults lie inside the 68% contour of the CIRS mass fimction. 
To compute the constraints in Figure 3 we assumed h=0.7 
and computed F from Equation (2). Figure 4 shows the 
effects of assuming different values of h or of assuming a 
fixed value of r=0.21. Clearly, the mass function com- 
puted with virial masses is now able to provide cosmolog- 
ical constraints competitive with those from X-ray mass 
estimates and even WMAP. 



3.2. Comparison With Other Methods 

We now compare the CIRS mass function and the result- 
ing constraints on cosmological parameters with previous 
studies. 

3.2.1. Mass Function 

Figure 5 shows the CIRS cluster mass function com- 
pared to other methods. HIFLUGCS uses the temperature 
and deprojected distribution of X-ray gas to estimate clus- 
ter masses from the ROSAT All-Sky Survey (Reiprich & 
Bohringer 2002). Bahcall et al. (2003) use a mass-richness 
relation to estimate the masses of clusters from the Early 
Data Release of SDSS. The CIRS mass function agrees 
well with both of these independent methods, indicating 
that the cluster mass function is robustly estimated at low 
redshift. 

3.2.2. Cosmological Constraints 

Figures 6-9 show a variety of other constraints in the 
(l^TOjCTs) plane. The colored contours in each Figure are 
for the CIRS virial mass function with the 68%, 95%, and 
99.7% confidence levels for flm and as- The two sets of 
contours shown in each panel show the 68% and 95% con- 
fidence levels from WMAP3. The dashed contours in Fig- 
ure 6 shows these levels for WMAPl, and the dash-dotted 
contours are the confidence levels from CFHTLS (Hoek- 
stra et al. 2005). Note that the contours for WMAPl are 
adapted from Figure 1 of Spergel et al. (2006), which dis- 
plays the contours in the flmh^-i^s plane. We increase 
the uncertainties in to account for the uncertainty 
in h. The 68% confidence levels of WMAP3 and CIRS 
overlap each other, demonstrating excellent agreement be- 
tween these two very different types of constraints. The 
CIRS contours lie closer to the WMAP3 contours than the 
CFHTLS results. 

Figure 7 compares the CIRS results with various lensing 
constraints. Red lines show the range of cosmic shear mea- 
surements, from the "low" result of Heymans et al. (2005, 
dash-triple-dotted line) to the "high" result of Massey 
et al. (2005, long-dashed line). The open square with ex- 
tremely small errorbars is the result of Seljak et al. (2006) 
for joint constraints from WMAP3, small-scale CMB mea- 
surements, SDSS galaxy clustering, supernovae, and the 
Lyman-a forest. The 68% confidence contour of the CIRS 
virial mass function encloses this point. Remarkably, the 
cosmic shear estimate from the CFHT Legacy Survey 
"Deep" sample (Semboloni et al. 2006, thick red short- 
dashed line in Figure 7) lies within the CIRS mass func- 
tion 68% contour, and the CFHTLS "Wide" sample over- 
laps significantly with the CIRS contours. The red dotted 
line shows the constraints from the mass function inferred 
from the M — Lx relation obtained from weak lensing 
measurements in SDSS (Viana et al. 2002). Finally, the 
thick red dash-dotted line shows the recent constraint of 
Dahle (2006) using weak lensing mass estimates of X-ray 
selected galaxy clusters. These constraints lie somewhat 
below the CIRS constraints and closer to the WMAP3. 
One potential difference between the samples is that the 
lensing clusters are at z=0.15-0.3, so cluster evolution may 
cause some of this difference. 

Blue lines in Figure 8 show constraints from X-ray mea- 
surements. The lower short-dashed line is the tempera- 



6 



ture function constraint of Seljak (2002); the upper short- 
dashed Une shows this relation increased by 20% in erg 
to show the systematic offset suggested by Rasia et al. 
(2005). The long-dashed line is from the X-ray mass func- 
tion of Rciprich & Bohringcr (2002), and the dash-dotted 
line shows the constraint from the X-ray luminosity func- 
tion (Allen et al. 2003). The diamond shows a recent X-ray 
measurements from Voevodkin & Vikhlinin (2004) using 
the baryon mass function and baryon fraction. The trian- 
gle shows the constraint from Picrpaoli et al. (2003) using 
the X-ray temperature function. The square at large flm is 
from combining the cluster abundance with the observed 
clustering (Schuecker et al. 2003, the two sets of errorbars 
show statistical and systematic errors). All error bars de- 
note 68% uncertainties. 

The upper dotted black line in Figure 9 is the virial 
mass function of Girardi et al. (1998b), and the thick 
black short-dashed line (and black cross) is from the mass- 
richness relation (Bahcall & Bode 2003). The thick triple- 
dot-dashed line is the constraint from Eke et al. (2005) 
for galaxy groups in the 2dFGRS. This line passes close 
to the middle of the CIRS 68% contour. The black long- 
dashed line shows the constraints of Tinker et al. (2005) 
from modeling the halo occupation distribution of galaxies 
and cluster mass-to-light ratios. Finally, the dash-dotted 
line shows the constraint from combining peculiar velocity 
surveys (Pike & Hudson 2005). 

The constraints on {Q,n,o'8) from the CIRS virial mass 
function are in excellent agreement with previous results. 
They lie between the upper and lower limits of cosmic 
shear estimates, and very close to the CFHTLS con- 
straints. Compared to optical estimates, CIRS lies slightly 
above Bahcall & Bode (2003), but agrees extremely well 
with Eke et al. (2005), who used a simplified virial-like 
equation to estimate masses. The CIRS constraints yield 
slightly larger values of {il„i,as) than most X-ray esti- 
mates, but the systematic temperature underestimate sug- 
gested by Rasia et al. (2005) would bring these results into 
reasonable agreement with the CIRS constraints. Finally, 
the constraint of Seljak et al. (2006) combines WMAP3 
with a number of independent probes to obtain estimates 
of flm and as with very small statistical uncertainties. The 
values of these parameters are shifted noticeably relative 
to the constraints from WMAP3 alone, and the sense of 
this shift is exactly in the direction of yielding better agree- 
ment with the CIRS mass function. 

3.3. Testing for Systematics 

It is tempting to conclude from the above comparisons 
that the CIRS virial mass function is an unbiased, reliable 
estimator of Sim and erg. Below, we overcome this tempta- 
tion and discuss potential systematic effects which could 
bias our results. 

3.3.1. Completeness of CIRS 

The CIRS sample is derived from four independent X- 
ray cluster catalogs constructed from the ROSAT All-Sky 
Survey. None of these surveys is complete: BCS/eBCS 
is estimated to be 80-90% complete, the main NORAS 
catalog is estimated to be 52% complete, the NORAS ex- 
tended sample in the 9''-14'' {5 >0) region is 82% com- 
plete (both NORAS estimates are obtained by comparing 



soiirce counts to the deeper REFLEX survey and assum- 
ing identical log N-log S distributions in the northern and 
southern hemispheres). XBACS searched for extended X- 
ray emission around all Abell clusters, so the CIRS sam- 
ple should include all Abell and ACO clusters (Abell et al. 
1989) with significant X-ray emission (although large off- 
sets between Abell and X-ray centers may cause XBACS 
to miss some systems). To estimate the approximate in- 
completeness in CIRS, we note that 53 of the 66 CIRS 
clusters arc in the northern hemisphere (5 >0) and thus 
in the area of sky covered by NORAS. Of these 53 clus- 
ters, 29 are NORAS (main catalog) clusters meeting our 
flux limit. If we assume that the 52% completeness of the 
main NORAS catalog applies to our redshift-limited sub- 
sample, there should be (29/0.52)=56 northern clusters 
in CIRS. Thus, the completeness of the composite CIRS 
catalog is «95%. 

We refit the mass function using only clusters in the 
REFLEX and NORAS main catalogs (omitting those in 
the NORAS 9^-14'' sample). The REFLEX southern cat- 
alog covers a small fraction of the DR4 area; including 
it eliminates the need to recalculate the effective area of 
the survey The REFLEX/NORAS sample should have a 
completeness similar to, but slightly larger than, the 52% 
completeness of the NORAS main sample. The cosmolog- 
ical constraints from this subsample assuming either 52% 
or 60% completeness are very similar to the constraints 
from the entire sample, although the confidence contours 
are larger. This result suggests that the use of the com- 
posite catalog of X-ray clusters introduces no significant 
bias in the mass fimction. 

Popesso et al. (2004) searched several X-ray catalogs 
of clusters and groups including the unpublished NO- 
RAS2 catalog, a deeper and more complete version of 
NORAS. They find four additional systems in the SDSS 
DR3 spectroscopic area which are not in any of the X- 
ray cluster catalogs used to construct the CIRS sample. 
Of these four systems, one (A1139) is included in the 
NORAS 9''-14'* catalog but hes below our flux hmit, an- 
other (NGC6173=A2197E) is sometimes included as part 
of A2197W, and Rines et al. (2002) find a fiux for A2197E 
that is below the CIRS flux limit. The remaining two sys- 
tems are apparently X-ray groups with velocity dispersions 
of (Tp<450km s"^ and virial masses <1 x lO^/i-^M©. This 
confirms that many of the systems missing from the CIRS 
sample are of relatively low mass and would probably lie 
below our adopted mass cutoff. 

We apply the V/Vmax test (Schmidt 1968) to the CIRS 
mass function sample to test for incompleteness. We find 
<V"/Knaa;>=0.518±0.035, Compared to the expected value 
of 0.5 for a uniform sample. This test shows that the CIRS 
sample is a representative sample. 

We recalculate the cluster mass function using a higher 
flux limit than the fiducial value of 3 x lO^^^crg s~^cm~^ 
(0.1-2.4 keV) cited for the parent X-ray catalogs. The 
mass function computed from the 46 clusters brighter than 
/x > 5 X 10~^^erg s~^cm~^ (0.1-2.4 keV) is almost iden- 
tical to our standard mass function. Figure 10 shows the 
effect of adopting this higher flux limit on our cosmological 
constraints (dash-dotted contours). The confldence con- 
tours get larger, but the best-fit parameters do not change 
significantly. 
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The minimum redshift we imposed to select the CIRS 
mass function sample reduces the number of low-mass 
systems sampled. Including these systems causes a sub- 
stantial increase in the number of systems at masses 

<lx IO^'^/i^^A/q, confirming our assertion that cosmic 
variance becomes significant in this mass range for the 
CIRS sample. 

Distinguishing a single, unrelaxed cluster from two 
merging clusters is a difficult observational problem. The 
CIRS sample contains several examples of multiple clus- 
ters (Rines & Diaferio 2006). We test the sensitivity of the 
CIRS mass function to this effect by making the extreme 
assumptions that either none or all multiple peaks are sep- 
arate halos. The affected systems are A1035A, A1035B, 
A1291A, A1291B, A1750A, A1750B, A2197W, A2197E, 
NGC5813, and NGC4636. Figure 10 shows that including 
all of these systems as well as the low-redshift systems (we 
call this the "maximal" CIRS mass function) makes little 
difference to our constraints in the (Qrm'^s) plane. 

A related completeness problem stems from the fact 
that, although we use galaxy dynamics to estimate masses, 
the CIRS sample is selected by X-ray flux (with an upper 
limit on redshift). This selection may cause the sample 
to underrepresent clusters with unusually low X-ray lumi- 
nosities for their masses. This incompleteness almost cer- 
tainly causes some of the deviations from the theoretical 
mass functions seen at masses below xlO^'^h~^MQ. How- 
ever, provided that the scatter in the Lx — M20Q relation is 
not too large, this incompleteness should be corrected for 
in the '^/Vmax{Lx) weighting. Also, our redshift cutoff of 
z < 0.1 likely increases the completeness because many of 
the lower-mass clusters are in the Abell catalog and some 
of the X-ray cluster catalogs explicitly searched for X-ray 
emission around Abell clusters. This effect likely domi- 
nates the competing effect that low-redshift, lower-mass 
clusters may have lower X-ray surface brightnesses. 

Even if the above effects all combine to bias our mass 
function to erroneously low number densities, we estimate 
a maximum incompleteness of a factor of two. Note that 
such a large incompleteness would yield best-fit parame- 
ters O„j=0.33 and (T8=0.87, outside the 99.7% confidence 
contour of our standard estimate. The best-fit contours 
computed assuming the NORAS completeness of 52% for 
the entire CIRS sample (i.e., approximately doubling the 
number densities) produces contours shifted towards larger 
values of Cl„i and, to a lesser degree, towards larger values 
of (78. Realistically, the incompleteness is probably larger 
at smaller masses. Correcting for such a mass-dependent 
incompleteness would tend to decrease the inferred esti- 
mate of as with a much smaller corresponding increase in 



3.3.2. Mass Estimation and Velocity Segregation 

The largest potential systematic effect we find is arti- 
ficially increasing the uncertainty in the virial mass esti- 
mates. As the assumed uncertainty increases, the con- 
volved mass function is more dramatically flattened rela- 
tive to the true mass function. Figure 11 shows the effects 
of increasing the virial mass uncertainty crj^g^ by a fac- 
tor of two (26% uncertainty, dashed contours) and by a 
factor of five (factor of two uncertainty, dash-dotted con- 
tours) without using the Vmax weighting of Equation 3. 



We find that artificially increasing the mass uncertainty 
by a factor of five and eliminating the Vmax weighting is 
sufficient to bring the CIRS mass function into agreement 
with the WMAP3 parameters, while the change is much 
smaller when the Vmax weighting is included. The effect 
of mass uncertainties is much more severe for the poorly 
sampled systems of Eke et al. (2005) and require detailed 
corrections from mock catalogs. Because the CIRS clus- 
ters are much more densely sampled, the effect should be 
miicli smaller for the CIRS mass function. Indeed, Figure 
9 shows that the constraints from the CIRS mass function 
agree well with those of Eke et al. (2005). 

We investigate two potential sources of systematic error 
caused by the mass range of our sample. Because the con- 
centration c is a weak function of mass (c oc M~°'^^), this 
may affect the comparison of the mass function via the 
conversion from Misob to M200. We refit the mass func- 
tion using the relation of Bullock et al. (2001) and find 
essentially the same results as for fixed c. Note that the 
CIRS sample showed little evidence of even a weak trend 
of c with mass (Rines & Diaferio 2006, but see Vikhlinin et 
al. 2006). Second, the uncertainty in the estimated mass 
may depend on the mass of the system. Less massive sys- 
tems are expected to contain fewer galaxies and thus may 
have larger statistical uncertainties. We find such a trend 
in the CIRS virial masses and redo the fits to the mass 
function assuming an uncertainty that declines with mass 
according to ay^^j^j = 0.34 — 0.021ogM in the convolution 

of the theoretical mass functions. The resulting confidence 
levels have nearly identical shapes and locations to our de- 
fault contours. 

The issue of velocity segregation between galaxies and 
dark matter remains an unresolved problem (Kauffmann 
et al. 1999a. b; Dicmand et al. 2004; Faltenbacher et al. 
2005). Velocity bias might be significant in the centers 
of clusters, where the galaxies have experienced dynami- 
cal friction. Dynamical friction should produce a smaller 
velocity dispersion for cluster galaxies than dark matter 
(Klypin et al. 1999; Yoshikawa et al. 2003; van den Bosch 
et al. 2005), although they may undergo more frequent 
mergers (due to their smaller velocities) or tidal disrup- 
tion, resulting in a larger velocity dispersion of the sur- 
viving cluster galaxies relative to the dark matter (Colin 
et al. 1999; Diemand et al. 2004; Faltenbacher et al. 2005). 
Note also that the velocity bias may depend on the mass of 
the cluster (Berlind et al. 2003) or the luminosities of the 
tracer galaxies (Slosar et al. 2006). An important caveat 
is that the subhalo distribution of cluster galaxies in simu- 
lations does not match the observed distributions of clus- 
ter galaxies: real cluster galaxies approximately trace the 
dark matter profile (e.g., Carlberg et al. 1997; Rines et al. 
2001; Biviano & Girardi 2003; Katgert et al. 2004; Lokas 
& Mamon 2003; Lin et al. 2004; Rines et al. 2004), while 
simulated cluster galaxies are significantly antibiased in 
cluster centers. This difference is likely caused by the re- 
silience of the stellar component of a subhalo (i.e., real 
galaxies) against tidal disruption relative to the total sub- 
halo in simulations, implying that a better understanding 
of the difference between galaxies and subhaloes is required 
for accurate models of velocity segregation (e.g., Diemand 
et al. 2004; Faltenbacher et al. 2006). For instance, the 
details of velocity bias depend on the selection used to 
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identify subhaloes, e.g., whether subhalo mass is measured 
prior to or after a subhalo enters the cluster (Faltenbacher 
et al. 2006). 

The amount of velocity bias in simulations is typically 
^10% for cluster-size haloes, although disagreement re- 
mains on whether this bias is positive or negative. Ve- 
locity bias of this size would lead to virial masses overesti- 
mated or underestimated by ~20%. Future work is clearly 
needed to determine the nature and significance of velocity 
bias for cluster mass estimates. There is no clear evidence 
from observations in favor of or against velocity bias in 
clusters, although the generally good agreement between 
X-ray, lensing, and virial mass estimates suggests that any 
velocity bias is not large (<20%; e.g., Girardi et al. 1998b; 
Popesso et al. 2004; Diaferio et al. 2005). Because sig- 
nificant velocity bias would produce incorrect virial mass 
estimates, the agreement between virial and caustic mass 
estimates is no guarantee that the caustic mass estimates 
are accurate. However, the mechanisms which cause ve- 
locity bias are less effective in the outskirts of clusters, so 
the caustic technique should be less afi^ected by velocity 
bias than estimates based on Jeans' analysis. 

Biviano ct al. (2006) test the robustness of virial mass es- 
timates in a cosmological hydrodynamic simulation. They 
find that virial mass estimates are generally reliable for 
densely sampled clusters, yielding overestimates of «10% 
for systems with more than 60 members. In general, they 
find that the most significant problem with c;stiniating 
virial masses is adequate removal of interlopers. The use 
of the caustic technique, a relatively conservative member- 
ship definition, should mitigate this problem for the CIRS 
sample. 

If one is willing to take the WMAP results as a prior, 
then one can estimate the sense and degree of velocity seg- 
regation by requiring that the CIRS mass function repro- 
duce the predictions of WMAP. The simplest implemen- 
tation of this method is to assume that the virial masses 
are biased by a (mass-independent) factor S such that 
the true cluster mass is given by M2oo,true = SM2oo,obs- 
If we use the WMAPl parameters as a prior, we ob- 
tain S = 1.14 ± 0.12, implying a velocity segregation of 
crgxy/(JDM «0.94±0.05. Taking the WMAP3 parameters 
as a prior, we obtain S = 0.61 ± 0.05, implying a veloc- 
ity segregation of agxy/crjjM «1.28±0.06. Seljak et al. 
(2006) combined constraints from WMAP3 and a variety 
of other tcchniriues (Lyman-a forest, SDSS galaxy cluster- 
ing, SDSS baryonic acoustic oscillations, supernovae, and 
smaller-scale CMB measurements) to improve constraints 
on cosmological parameters. They find that the best-fit 
estimate of as increases to as = 0.847 ± 0.022 with essen- 
tially no change in = 0.235 ± 0.008. If we repeat the 
above analysis using these parameters as priors, we find 
S = 0.83 ± 0.08 and ag^y/aoM wl.l0±0.05. Thus, an 
independent measurement of cosmological parameters can 
constrain the detailed properties of cluster galaxies (Eke 
et al. 2005; Berhnd et al. 2006). 

Blue and red galaxies in clusters have dramatically dif- 
ferent spatial distributions. By using all galaxies within 
the caustics to estimate virial masses, we may overesti- 
mate M200- We recalculate the virial mass profiles using 
only red galaxies, where we define red galaxies as those 
with (5-r)>-0.05(Mr-h20)-|-0.7 where Mr is the absolute 



magnitude in r band and g and r are Petrosian magni- 
tudes. We find that there is very little difference in the re- 
sulting virial mass estimates and, as a result, little change 
in the estimated mass function or cosmological parameters 
(Figure 3). This result differs from the predictions of Bi- 
viano et al. (2006) that the virial masses estimated with 
all galaxies are biased high by a factor of ~1.40, while the 
virial masses computed with red galaxies only are approx- 
imately unbiased {Mmr,red ~ l-01Ait,ir,trMe)- We suggest 
that this difference is due to the difference in the contrast 
of caustics in observations and in simulations (Rines et al. 
2003; Rines & Diaferio 2006). A more detailed study of 
the velocity distribution of red and blue galaxies in and 
around clusters suggests that intrinsic differences in the ve- 
locity dispersions of the two populations are substantially 
reduced when using the caustic technique to select mem- 
bers (Rines & Diaferio in prep.). This study also shows 
that, when using the caustic technique, the inclusion or 
exclusion of blue galaxies produces little change in the ve- 
locity dispersion of the total population (in part because 
blue galaxies are rare in clusters). 

Most of the CIRS clusters are sampled only once by 
SDSS. Thus, the fiber placement constraints of SDSS (no 
two fibers can be placed closer than 60") causes non- 
uniform sampling of the cluster galaxies. Biviano et al. 
(2006) show that non-uniform sampling of cluster galax- 
ies by multislit or multifiber surveys may bias virial mass 
estimates (the cluster core is imdersampled so the virial 
radius is biased). Unfortunately, few clusters have suf- 
ficient coverage in both SDSS and elsewhere to test this 
effect directly. Followup spectroscopy of the CIRS clusters 
to complete the spectroscopic samples from SDSS would 
test the significance of this effect. 

3.3.3. Peculiar Velocities 

We make no corrections for peculiar velocities of the 
CIRS clusters because no reliable peculiar velocity esti- 
mates exist for most of the sample. The closest clusters 
will be most strongly affected by this approach, because 
their distance estimates are most sensitive to a given cor- 
rection for peculiar velocity. Our minimum rcdshift of 
z>0.02 therefore significantly reduces this uncertainty. 

We estimate the magnitude of this effect by adding a 
random peculiar velocity from a Gaussian distribution of 
width 400km s"^ (Sheth fc Diaferio 2001) to each cluster 
rcdshift and recomputing the cluster mass function and 
cosmological constraints. Several trials using this proce- 
dure find confidence contours for {Qrm'^s) that agree well 
with our default contours. We conclude that peculiar ve- 
locities have a small impact on our results, although ex- 
tensive Monte Carlo simulations could better quantify the 
resulting systematic uncertainty. 

3.3.4. Large-Scale Structure 

Large-scale structure could have a significant impact on 
our mass function if the scale of overdcnsities is large com- 
pared to the volume of the survey. The CIRS sky coverage 
is that of the spectroscopic survey of SDSS DR4. This 
means that CIRS does not include the Coma cluster or 
much of the Great WaU (Geller & Huchra 1989), but it 
does include the Sloan Great Wall (Gott et al. 2005). Be- 
cause of the redshift limit, the maximum volume probed is 

10^ h~^Mpc^ (Figure 1). Figure 12 shows a histogram of 
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the redshift ditribution of the CIRS mass funetion sample. 
The sohd curve in Figure 12 shows the expected redshift 
distribution for the CIRS sample given the X-ray luminos- 
ity ftmction of Bohringer et al. (2002). There is no strong 
evidence for large-scale variations in this histogram. It is 
difficult to determine whether the CIRS sample contains 
an overdense or underdense sample of the universe. The 
eventual release of SDSS DR5 and the SDSS-II Legacy 
Survey will approximately double the spectroscopic area 
covered, so those data will provide a good test of the rep- 
resentativeness of the CIRS sample. 

3.3.5. Scaling Relations 

An alternative method of estimating mass functions is 
to measure simple observables such as X-ray luminosity 
or galaxy richness and then apply scaling relations to es- 
timate cluster masses. Although this method has been 
applied extensively (Picrpaoli et al. 2001; Seljak 2002; 
Pierpaoli et al. 2003), the reliability of the mass function 
and/or cosmological constraints is limited by the reliabil- 
ity of the scaling relations. Scaling relations often have 
large scatter (Popesso et al. 2004) and they are difficult 
to caUbrate (Evrard et al. 1996; Stanek et al. 2006). We 
recalculate the CIRS mass function using the observed X- 
ray luminosities and two scaling relations for Lx — Af200j 
one based on virial masses (Rines & Diaferio 2006) and one 
based on masses from X-ray temperatures (Popesso et al. 
2004). The resulting mass function for the latter relation 
contains many fewer systems at M<10^'^h~^ Mq than the 
CIRS mass function (Figure 5). This difference seems to 
be caused by a difference in the scaling relations deter- 
mined by Popesso et al. (2004) and Reiprich & Bohringer 
(2002). Rines & Diaferio (2006) noted that their Myir-Lx 
relation was offset from that of Popesso et al. (2004), al- 
though the M — Tx relations agree. The Lx — M200 of 
Reiprich & Bohringer (2002) and Rines & Diaferio (2006) 
agree with each other, and so do the resulting estimates of 
the mass function. We therefore caution that great care 
and judgment must be taken when using scaling relations 
to measure the cluster mass function. 

3.3.6. Sample Selection 

There are several methods of identifying clusters. Four 
techniques are used to identify clusters: galaxy overden- 
sities (either in projection or in redshift space), identi- 
fying extended X-ray sources, identifying peaks in gravi- 
tational lensing maps, and identifying Sunyaev-Zeldovich 
decrements in maps of the cosmic microwave background. 
Lensing and SZ surveys are both very promising but have 
yet to produce large, well-selected samples. 

Identification of galaxy overdensities has proved very 
successful, but this method is subject to projection ef- 
fects, especially when performed without redshift infor- 
mation. Eke et al. (2005) study the mass function from an 
optically-selected group catalog containing systems with 
as few as two members. They justify this procedure by 
extensive comparison to mock catalogs, but the resulting 
cosmological constraints then have additional systematic 
uncertainty due to any differences between the mock cat- 
alogs and the observations. 

Popesso et al. (2006) describe a targeted search for X- 
ray emission around Abell clusters in SDSS. They find that 



many Abell clusters with apparently Gaussian velocity dis- 
tributions lie far below the Lx — M relation of X-ray se- 
lected clusters, although they lie on the same L^pt—M rela- 
tion. If these X-ray faint clusters are indeed very massive, 
then the constraints on Q,m and erg from the CIRS mass 
function are underestimates. Note, however, that these 
changes would lead to larger disagreements with other esti- 
mates (Figure 6), unless either velocity segregation is large 
(and therefore the masses are overestimated) or the true 
uncertainty in the measured virial masses is much larger 
than the statistical uncertainty (or some combination of 
these effects). 

It would be very instructive to compare the CIRS mass 
function to a cluster mass function constructed with op- 
tical cluster catalogs from SDSS (e.g., Miller et al. 2005), 
but this analysis lies outside the scope of this work. De- 
tailed comparisons of cluster catalogs selected with differ- 
ent techniques (e.g., Donahue et al. 2002) are important 
for using clusters to constrain cosmological parameters. 

3.3.7. Uncertainties in the Predicted Mass Function 

Most early attempts to measure cosmological param- 
eters from the observed mass function of clusters used 
Press-Schechter formalism (Press & Schechter 1974). In 
recent years, however, numerical simulations have achieved 
sufficient resolution to directly measure the mass function 
of dark matter haloes in cosmologically representative vol- 
umes. Sheth & Tormen (1999) proposed a modified form of 
the Press-Schechter relation that provides a better match 
to simulations. Jenkins et al. (2001) provide a third form of 
this relation which provides a good match to simulations of 
the Virgo consortium. We use their fitting formula for the 
mass function along with the power spectrum of Bardeen 
et al. (1986). Reiprich & Bohringer (2002) notes that the 
best-fit cosmological parameters are unchanged if one uses 
the power spectrum of Eisenstein & Hu (1999) instead. 
Warren et al. (2005) use 16 nested simulations to test the 
effects of numerical resolution on the mass function. Al- 
though they find significant differences from the Jenkins 
et al. (2001) relation for extreme masses {\i<\Q^^h~^MQ 
andM> 2x 10^^/i"^Mq), the Jenkins et al. (2001) relation 
provides a good approximation over the mass range consid- 
ered here. The constraints in the {Q.m,<y») plane obtained 
by using the Press-Schechter or Sheth- Tormen form are 
nearly identical to those obtained with the Jenkins et al. 
(2001) relation. 

White (2002) discusses the importance of carefully defin- 
ing the mass of a halo when comparing theoretical and 
observed mass functions. In particular, the universality 
of the mass function depends on the algorithm used to 
measure the mass. White (2002) concludes that the re- 
lation of Equation (9) of Jenkins et al. (2001) provides a 
good match to the mass measured within a spherical over- 
density of 180 times the background (not critical) density. 
Both White (2002) and Hu & Kravtsov (2003) recommend 
computing the mass function for a given set of cosmologi- 
cal parameters using the 180b definition of mass and then 
using an NFW profile with concentration C2oo=5 to con- 
vert the observed and theoretical mass scales to a common 
scale. Here we convert the mass scale of the theoretical 
mass function so that the models can be shown against 
the observations more easily. There appears to be con- 
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sensus in the literature that the systematic uncertainties 
in the predicted mass function for different cosmological 
models are of the order 15% (White 2002; Evrard et al. 
2002; Hu & Kravtsov 2003; Stanek et al. 2006). White 
(2002) suggests that higher accuracy may require a suite 
of cosmological simulations run with different values of 
and (78 with subsequent analysis using a variety of algo- 
rithms for defining the mass of a halo. 

Some authors have noted that the most massive clusters 
mcasiirc fluctuations on a scale larger than 8/?,^^Mpc. In 
particular, Pierpaoli et al. (2001) note that clusters with 
temperatures of 6.5 keV sample fluctuations on a scale of 
KlAh~^Mpc. The CIRS sample probes lower mass sys- 
tems than many previous studies: the median mass of 
the sample is 2xl0^'^h~^ Mq, corresponding to a scale of 
8.3/i~^Mpc for 17^=0.3. Thus, we constrain ag directly. 
Virial masses have an advantage over X-ray masses in this 
respect because the properties of the ICM for clusters with 
M^lO^^ft,~^M0 may be strongly affected by feedback from 
supernovae and AGN (e.g., Loewenstein 2000), and the 
ICM in groups may not be in hydrostatic equilibrium. 

4. VELOCITY DISPERSION FUNCTION 

The abundance of systems as a function of velocity dis- 
persion is an important input to models of gravitational 
lensing. Sheth et al. (2003) and Mitchell et al. (2005) com- 
pute the velocity dispersion function of early-type galaxies 
in SDSS. One might expect a smooth transition between 
the velocity dispersion function of galaxies and the veloc- 
ity dispersion function of groups and clusters (Shimasaku 
1993). Zabludoff et al. (1993) found the cluster velocity 
dispersion function for systems identifled in the CfA Red- 
shift Survey (see also Pisani et al. 2003). They found good 
agreement with the velocity dispersion function implied by 
the X-ray temperature function. Figure 13 compares the 
velocity dispersion functions of these two types of systems. 
The cluster velocity dispersion function shows a dcflcit of 
low-<j systems relative to the number needed to produce 
a smooth function. The dashed line in Figure 13 shows 
the velocity dispersion function for our "maximal" mass 
function sample including cluster pairs and low-redshift 
systems. This estimate is more; consistent with a smooth 
transition between dark matter haloes containing one ver- 
sus many galaxies, although the two types of systems still 
do not join smoothly. A careful study of the abundance 
of galaxy groups that span the intermediate mass and ve- 
locity dispersion range would shed light on this transition. 
Earlier studies (Zabludoff et al. 1993; Mazure et al. 1996; 
Girardi & Giuricin 2000; Heinamaki et al. 2003; Pisani 
et al. 2003) were generally limited to small, nearby vol- 
umes, so cosmic variance will continue to be an important 
limitation to these efforts. The much larger samples avail- 
able with 2dFGRS and SDSS should alleviate this problem 
(Martinez et al. 2002; Eke et al. 2005). 

Figure 14 shows the cumulative velocity dispersion func- 
tion of the CIRS sample. Our estimate is lower than that 
of Pisani et al. (2003) from the CfA2 Redshift Survey, but 
they note that their estimate is larger than most previ- 
ous estimates, a result they attribute to cosmic variance 
due to the presence of the Great Wall in the survey. The 
CIRS velocity dispersion function agrees with these earlier 
estimates (e.g., Zabludoff et al. 1993; Mazure et al. 1996). 



However, Casagrande & Diaferio (2006) analyze the prop- 
erties of groups in the Updated Zwicky Catalog (UZC) and 
flnd similar results to Pisani et al. (2003). 

5. DISCUSSION AND FUTURE WORK 

We use the Fourth Data Release of the Sloan Digital 
Sky Survey to measure the mass function of galaxy clus- 
ters in the local universe. We select clusters based on their 
X-ray flux and a redshift limit (z < 0.1). X-ray selection 
has the crucial advantage over optical selection of yield- 
ing a well-defined selection function. We measure cluster 
masses from the dynamics of their member galaxies using 
both the virial theorem and the caustic technique. The 
resulting mass functions agree well with previous deter- 
minations using X-ray properties and the mass-richness 
relation. We c;onclude that the cluster mass function is 
measured robustly in the local universe. We compare this 
mass function to predictions of various cosmological mod- 
els to obtain constraints in the {il„i,(7s) plane. With the 
CIRS virial mass function (and assuming h=0.7), we es- 
timate flrri=0.24to[lt and a-8=0.92t^;?^, or <78=0.84±0.03 
when holding 1^^=0.3 fixed. 

The CIRS virial mass function agrees well with the 
WMAPl and lies slightly above the WMAP3 results, al- 
though still within the 95% confidence levels. This result 
contrasts slightly with Reiprich (2006), who finds that the 
HIFLUGCS mass function lies much closer to WMAP3 
than WMAPl. However, other investigators have found 
mass functions using X-ray data that lie closer to the CIRS 
mass function, and there remain open questions about the 
systematic uncertainties in the X-ray mass estimates (Ra- 
sia et al. 2005; Stanek et al. 2006). Virial masses have the 
following advantages: they measure the masses at larger 
radii, where less extrapolation is required to compare with 
theoretical results, they are less sensitive to the complex 
physics in the centers of clusters, and they can be applied 
to lower mass systems for which nongravitational physics 
greatly complicates X-ray mass estimates, thus eliminat- 
ing any uncertainty due to possible scale dependence of 
the estimate of erg (the clusters of, e.g., Pierpaoli et al. 
2001, sample fluctuations on larger scales, an, which must 
then be converted to erg). We investigate several potential 
sources of systematic error and conclude that the most sig- 
nificant effect would be a large underestimate of the uncer- 
tainties in the virial mass estimates. Recent investigations 
with numerical simulations Biviano et al. (2006) suggest 
that the scatter in virial mass estimates of well-sampled 
clusters is smaller than the amount required to bring our 
results into agreement with the WMAP3 parameters, al- 
though further confirmation of this result would be useful. 

Although obtaining virial mass estimates from galaxy 
dynamics for large numbers of clusters requires much ob- 
servational effort, the reward is an independent measure- 
ment of the mass function and a powerful confirmation of 
the current cosmological model. This work further demon- 
strates the power of the mass function as a fundamen- 
tal constraint on cosmological parameters. The SDSS-II 
Legacy Survey should substantially enlarge the sample of 
local clusters and reduce both statistical and systematic 
uncertainties in the mass function. 

The CIRS mass function provides a critical normaliza- 
tion of the mass function in the nearby universe. An accu- 
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rate normalization of the local mass function is a prereq- 
uisite for future studies of the evolution of the mass func- 
tion. Future determinations of the mass function at higher 
redshifts will directly probe the growth rate of structure 
and provide further constraints on cosmological parame- 
ters (e.g., Haiman et al. 2001; Hu 2003; Vikhlinin et al. 
2003; Henry 2004; Majumdar & Mohr 2004; Voevodkin & 
Vikhlinin 2004). We conclude that a large spectroscopic 
program to measure virial masses of an X-ray selected sam- 
ple at inodcratc rcdshift would provide an independent 
measurement of the evolution of the mass function. 
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Fig. 1. — The maximum volume sampled by CIRS as a function of X-ray luminosity. The top axis shows the equivalent mass assuming the 
Lx — M200 relation of the CIRS sample. 
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Fig. 2. — The mass function of the CIRS sample. Red and blue lines are computed using virial masses and caustic masses respectively. The 
thick dash-dotted lines show the mass functions computed using the cosmological parameters from the WMAPl results (upper) and WMAP3 
results (lower) using the results of Jenkins et al. (2001). The red dashed line shows the best-fit mass function for the CIRS virial mass 
function. The vertical line indicates the minimum mass we use to constrain cosmological parameters. The dotted lines and errorbars show 
the mass function computed after removing the minimum rcdshift and including all possible mergers as separate systems. This demonstrates 
the importance of cosmic variance at these low masses. 
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Fig. 5. — Cumulative mass function of the CIRS sample. Red and blue lines are computed using virial masses and caustic masses respectively. 
Open circles show the HIFLUGCS mass function computed with X-ray masses (Reiprich & Bohringer 2002) and triangles show the mass 
function computed from early SDSS data using a mass-richness relation (Bahcall et al. 2003). Green and magenta lines show the mass function 



of the CIRS sample computed with the Lx 
respectively. 



Mx relation of Popesso et al. (2004) and the Lx — Myir relation of Rines & Diaferio (2006) 
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Fig. 6. — Cosmological constraints from the CIRS virial mass function compared to other results. Colored contours show l-2-3a confidence 
levels for Qrn and eg. The dashed, solid, and dash-dotted contours show the 68% and 95% confidence levels from WMAPl, WMAP3, and 
the CFHTLS Wide survey respectively Spergel et al. (adapted from Figures 1 and 7 of 2006). 
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Fig. 7. — Cosmological constraints from the CIRS virial mass function compared to other results. Colored contours show 1-2-3(j confidence 
levels for Qrn and erg. Dash-triple-dotted (Heymans et al. 2005) and long-dashed (Massey et al. 2005) lines show the range of constraints from 
cosmic shear measurements. The thick red short-dashed line shows the cosmic shear estimate from the CFHT Legacy Survey "Deep" sample 
(Semboloni et al. 2006). The red dotted line shows the constraints from the mass function inferred from the M — Lx relation obtained from 
weak lensing measurements in SDSS (Viana et al. 2002). The thick red dash-dotted line shows the constraint of Dahle (2006) using weak 
lensing mass estimates of X-ray selected galaxy clusters. The open square with extremely small errorbars is the result of Seljak et al. (2006) 
for joint constraints from WMAP3, small-scale CMB measurements, SDSS galaxy clustering, supernovae, and the Lyman-a forest. The solid 
contours extending to the lower left show the 68% and 95% confidence levels from WMAP3. 
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Fig. 8. — Cosmological constraints from the CIRS virial mass function compared to other results. Colored contours show 1-2-3(j confidence 
levels for f2m and erg . Blue lines show constraints from X-ray cluster data. The lower short-dashed line is the temperature function constraint 
of Seljak (2002); the upper short-dashed line shows this relation increased by 20% in trg to show the systematic offset suggested by Rasia 
et al. (2005). The long-dashed line is from the X-ray mass function of Reiprich & Bohringer (2002), and the dash-dotted line shows the 
constraint from the X-ray luminosity function (Allen et al. 2003). The diamond shows a recent X-ray measurements from Voevodkin & 
Vikhlinin (2004) using the baryon mass function and baryon fraction. The triangle shows the constraint from Pierpaoli et al. (2003) using the 
X-ray temperature function. The square at large Qrn, is from combining the cluster abundance with the observed clustering (Schuecker et al. 
2003, the two sets of errorbars show statistical and systematic errors). All error bars denote 68% uncertainties. The solid contours extending 
to the lower left show the 68% and 95% confidence levels from WMAP3. 
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Fig. 9. — Cosmological constraints from the CIRS virial mass function compared to other results. Colored contours show 1-2-3(7 confidence 
levels for Qm, and eg. Black lines show constraints from optical cluster data. The upper dotted black line is the virial mass function of 
Girardi et al. (1998b), and the thick black short-dashed line (and black cross) is from the mass-richness relation (Bahcall & Bode 2003). The 
triple-dot-dashed line is the constraint from Eke et al. (2005) for galaxy groups in the 2dFGRS. The long-dashed line shows the constraints 
of Tinker et al. (2005) from modeling the halo occupation distribution of galaxies and cluster mass-to-light ratios. The red dash-dotted line 
shows the constraint from combining peculiar velocity surveys (Pike & Hudson 2005). All error bars denote 68% uncertainties. The solid 
contours extending to the lower left show the 68% and 95% confidence levels from WMAP3. 
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Fig. 11. — Cosmological constraints from the CIRS virial mass function. Solid contours show l-3a confidence levels for Qm and crs- Dashed 
and dash-dotted contours show the l-3o" confidence levels from adopting a mean uncertainty 2x and 5x the statistical uncertainty respectively. 
The solid contours extending to the lower left show the 68% and 95% confidence levels from WMAP3. 
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Fig. 12. — Rcdshift histogram of the CIRS mass function sample (thick black line). The thin black histogram shows the redshift histogram 
of the "maximal" CIRS sample including cluster pairs and systems at low redshift. The solid curve shows the expected number of clusters 
for the X-ray luminosity function of Bohringer et al. (2002). 
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Fic;. 13. The velocity dispersion function of the CIRS sample. The solid curve shows the best-fit velocity dispersion function of early-type 
galaxies in SDSS and the black dots indicate the deviations from this fit at log(T>2.4 (Mitchell et al. 2005). The dash-dotted line is the 
velocity dispersion function of the "maximal" CIRS sample including cluster pairs and low-redshift systems. 
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Fig. 14. — The cumulative velocity dispersion function of the CIRS sample (squares). The triangles show the velocity dispersion function 
of the "maximal" CIRS sample including cluster pairs and low-redshift systems. The dash-dotted lines show the estimate and 95% confidence 
range from Pisani et al. (2003), and circles show their compilation of previous results at cr=750km s~^. Dashed and dotted lines show 
estimates from Zabludoff et al. (1993) and Mazure et al. (1996). 



